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1.  Introduction 

A.  Background 

The  development  of  practical  quantum  key  distribution  (QKD)  systems  depends  upon  the 
availability  of  efficient,  low-noise  single -photon  detectors.  An  important  class  of  QKD  systems 
is  designed  to  use  standard  telecommunications  optical  fiber  as  the  link  conveying  quantum 
information,  and  therefore  operation  at  a  wavelength  of  1310  or  1550  nm  is  required.  In  order 
to  detect  single  photons  at  these  wavelengths,  commercially- available  Ge  or  InGaAs  avalanche 
photodiodes  (APDs)  have  been  used,  albeit  in  a  manner  for  which  they  were  not  designed — gen¬ 
erally,  they  must  be  cooled  to  low  temperatures  (~77-150K)  and  operated  above  their  reverse 
breakdown  voltages  (“Geiger  mode”)  in  order  to  detect  single  photons  with  acceptably  high 
detection  efficiencies  and  acceptably  low  dark  count  rates. 

Under  the  extreme  operating  conditions  required  for  single-photon  detection,  there  is  a  wide 
variability  in  the  performance  of  detectors  from  manufacturer-to-manufacturer,  and  even  from 
device-to-device  made  by  the  same  manufacturer.  Moreover,  as  manufacturers  have  refined  suc¬ 
cessive  generations  of  APDs  in  order  to  improve  their  performance  in  standard  telecom  applica¬ 
tions,  they  have  actually  made  them  worse  for  single -photon  detection. 

A  further  very  severe  problem  is  that  specific  APDs  that  have  been  used  in  past  QKD  experi¬ 
ments  are  now  becoming  unavailable  to  researchers.  In  some  cases,  large -volume  orders  by 
optical  network  developers  are  consuming  the  entire  supply  of  some  devices;  others  are  no  longer 
available  as  discrete  detectors,  but  are  being  packaged  with  other  components  (e.g.,  GaAs  pre¬ 
amplifiers)  in  a  way  that  makes  them  unsuitable  for  use  as  single -photon  counters;  still  others  are 
being  discontinued  as  demand  declines  due  to  the  slowdown  in  the  telecom  industry. 

Hence,  in  order  to  continue  the  advancement  of  QKD  systems  and  to  enable  new  experi¬ 
ments  involving  transport  of  quantum  information  over  long  distances  via  optical  fiber,  it  is  essen¬ 
tial  to  develop  new  detectors  with  high  detection  efficiency  and  low  dark  count  rate. 

I  20050901  060 


1 


B.  Program  Objectives: 

1.  To  understand,  through  modeling  and  experimentation,  the  physical  parameters  that  are 

most  important  for  achieving  high-efficiency,  low-noise  detection  of  single  1310/1550  nm 
photons  using  avalanche  photodiodes. 

2.  Using  this  understanding,  to  design,  fabricate,  and  characterize  APDs  for  single-photon 

detection  at  telecom  wavelengths 

3.  To  test  these  improved  APDs  in  a  quantum  cryptography  system  and  understand  how  the 

improved  detector  characteristics  improve  the  performance  of  the  QKD  system. 

4.  To  investigate  nonlinear  frequency  conversion  combined  with  Si  APDs  as  an  alternative  to 

the  use  of  InGaAs  APDs  for  detection  of  single  telecom-wavelength  photons. 

C.  Program  Participants: 

IBM  Research 
Almaden  Research  Center 
Dr.  William  E  Risk  (PI) 

Dr.  Donald  S.  Bethune 

University  of  Texas  at  Austin 
Prof.  Joe  C.  Campbell  &.  graduate  students 

University  of  California  at  San  Diego 
Prof  Yu-Hwa  Lo  &  graduate  students 

D.  Program  Approaches: 

Under  this  program,  we  are  pursuing  three  independent  approaches  to  single -photon 
detection  at  telecom  wavelength,  which  are  summarized  here: 

1 .  Optimization  of  lowmoise  InGaAs  APDs  for  detection  of  single  photons  at  telecom  wave  - 
lengths 

This  part  of  the  program  aims  to  extend  earlier  work  done  by  Dr.  Joe  Campbell’s 
group  at  UT  Austin  on  ultra-low  noise  APDs  from  the  regime  of  conventional  telecom 
applications  into  the  single-photon  regime.  In  earlier  work,  the  UT  Austin  group  used 
APD  devices  with  very  thin  multiplication  regions  to  obtain  the  lowest  excess  noise  and 
the  highest  gain-bandwidth  products  reported  for  any  APD.  They  showed,  both  experi¬ 
mentally  [1-4]  and  theoretically  [5],  that  the  low  noise  is  due  to  the  non-local  nature  of 
the  impact  ionization  process  and  the  fact  that  non-locality  becomes  a  dominant  effect  in 
very  thin  high-field  regions.  In  addition,  they  achieved  ultra-low  noise  operation  by  using 
heterojunctions  to  provide  greater  localization  of  impact  ionization  than  can  be  achieved 
in  spatially  uniform  structures.  This  approach  they  call  “impact  ionization  engineering”  or 
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PE.  The  initial  demonstration  of  this  concept  used  the  GaAs/Al^Gaj^As  material  system; 

under  this  program,  the  concept 
is  being  extended  to  the  “long- 
wavelength”  materials  that  are 
widely  used  for  fiber  optic  com¬ 
munications. 


2.  Development  of  Si/InGaAs 
“Separate  Hetero  Interface 
Photodetectors" 

The  second  approach  to  high- 
efficiency,  low-noise  detection 
being  pursued  under  this  project 
is  the  “separate  hetero-interface 
photodetector  (SHIP) ,”  which 
combines  a  long-wavelength 
absorption  region  with  a  separate 
low-noise  multiplication  region  made  from  a  heterogeneous  material  (Figure  1).  For 
example,  this  type  of  structure  has  been  successfully  demonstrated  using  an  In^j^Ga^^^As 
absorber  and  a  Si  multiplication  region  joined  together  by  fusion  bonding  [6] .  We  believe 
that  devices  based  on  this  approach  may  be  capable  of  operating  as  single-photon  detec¬ 
tors  at  temperatures  near  room  temperature  or  requiring  only  modest  thermoelectric 
cooling. 

Earlier  work  on  InGaAs/Si  APDs  in  Dr  Yu-Hwa  Lo’s  group  at  UCSD  and  at  Nova- 
Crystals,  Inc.  (the  industrial  partner  of  the  UCSD  group)  has  focused  on  two  important 
aspects  for  optical  communications,  namely  high  gain-bandwidth  product  and  low  opera¬ 
tion  voltage.  Under  this  program,  the  emphasis  is  on  reduction  of  the  dark  current  and 
improvement  of  the  noise  figure.  It  was  expected  that  InGaAs/Si  heterostructure  APDs 
can  achieve  the  theoretical  limit  of  the  dark  current  with  (1)  improved  device  processing, 
(2)  improved  wafer  fusion  interface,  and  (3)  optimized  epitaxial  layer  structure.  Specific 
goals  of  this  part  of  the  project  at  the  outset  of  the  program  included  replacing  the  mesa 
structure  of  earlier  designs  with  a  planar  structure  with  guard  rings;  applying  a  low  tem¬ 
perature  fusion  process  in  high  vacuum  will  to  replace  the  previously- used  high  tempera¬ 
ture,  atmosphere  fusion  process;  and  using  an  optimized  doping  profile  to  minimize  the 
tunneling  current  and  to  facilitate  carrier  transport. 

3.  Nonlinear  frequency  upconversion  and  silicon  APDs 

A  third  approach  that  is  being  explored  for  detection  of  single  1.31/1.55-|a.m  photons 
at  or  near  room  temperature  uses  nonlinear  optical  techniques  to  upconvert  telecom 
wavelength  photons  to  a  visible  wavelength;  then  low- dark- count,  high  quantum  effi- 
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Figure  I.  Schematic  cross-section  of  a  separate  hetero-interface  photo¬ 
detector  (SHIP) 
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ciency  Si  APDs  can  be  used  to  detect  the  upconverted  photon  (Figure  2).  For  example, 
by  using  sum-frequency  generation  in  a  nonlinear  crystal,  it  is  possible  to  mix  a  single 
1.55-jj,m  photon  with  a  strong  1.3'|lm  pump  photon  to  create  a  single  photon  at  0.71  |a.m, 
a  wavelength  that  can  be  detected  efficiently  using  a  Si  APD.  Since  Si  APDs  are  (rela¬ 
tively)  insensitive  to  1.3/1.55  p,m,  it  should  be  possible  to  provide  sufficient  filtering  so 
that  only  the  generated  0.71-|lm  photon  triggers  an  avalanche. 

In  this  portion  of  the  program,  we  are  investigating  the  feasibility  of  this  approach 
and,  if  it  is  deemed  feasible,  we  will  work  toward  development  of  a  practical  prototype 
detection  system.  Planned  work  includes:  studying  the  inherent  sensitivity  of  Si  single - 
photon  counting  APDs  to  sub-bandgap  light  (A,<1  jlm)  to  evaluate  filtering  requirements. 


1310  nm 
single  pliolon 


Nonlinear  Filter 


single  photon 


Si  APD 


Figure  2.  Diagram  showing  the  concept  for  detect 
tion  of  single  telecom-wavelength  photons  using 
nonlinear  frequency  conversion. 
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evaluating  possible  upconversion  schemes  compatible  with  practical  nonlinear  materials 
(including  periodically-poled  materials  and  waveguide  geometries)  and  pump  sources,  and 
comparing  performance  of  such  a  system  to  commercially- available  Ge  and  InGaAs  APDs 
and  other  APDs  developed  under  this  program  as  described  above. 

Another  key  aspect  of  this  program  involves  characterization  of  the  detectors  that  are 
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developed  using  the  approaches  described  above.  Quantum  key  distribution  places  particular 
and  rather  specialized  demands  on  the  detectors;  “good”  performance  is  defined  rather  dif¬ 
ferently  than  in  conventional  telecom  applications.  It  is  therefore  necessary  to  character¬ 
ize  detectors  in  a  way  that  fairly  evaluates  the  characteristics  that  are  relevant  to  quantum 
cryptography.  Hence,  one  role  of  IBM  group  is  to  characterize  any  candidate  detectors  in  a 
manner  consistent  with  similar  measurements  that  they  previously  performed  on  commer¬ 
cially- available  detectors  and  to  evaluate  the  performance  of  promising  devices  in  the  IBM 
prototype  quantum  cryptography  system  (Figure  3)  [7]. 


IL  Third  Year  Progress 

A.  IBM  Almaden  Research  Center  Research  Prototype  Single  Photon  Detector 

1.  Introduction 

Detection  of  single  photons  at  telecom  wavelengths  for  quantum  cryptography  requires 
not  only  the  detector  element  itself  (typically  an  InGaAs  avalanche  photodiode),  but  also  a 
constellation  of  specialized  electronic  elements.  This  auxiliary  electroruc  circuitry  performs 
a  number  of  crucial  functions  (Figure  4).  It  generates  a  short  (~'l  ns)  bias  pulse  which, 
when  combined  with  a  DC  bias  voltage  applied  to  the  APD,  momentarily  biases  the  APD 
above  breakdown,  enabling  it  to  produce  an  avalanche  in  response  to  absorption  of  a  single 


Figure  4-  Bbck  diagram  showing  the  functions  that  must  be  performed  by  the  auxiliary  electronics  in  order 
to  enable  an  avalanche  photodiode  to  detect  single  photons. 


photon.  However,  applying  such  a  short  pulse  to  the  capacitive  diode  generates  a  transient 
current  that  can  exceed  the  current  resulting  from  the  avalanche.  Thus,  the  supporting 
electronics  must  also  separate  the  photon-induced  signal  from  the  capacitive  transient.  In 
the  IBM  prototype  QKD  system,  we  use  a  patented  balanced  cancellation  scheme  (US 
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Figure  5.  Laboratory 
apparatus  used  in  earlier 
detector  characterization 
work,  (a)  Electronic  gating 
implemented  with  Minicir- 
cuits  modules,  (b)  liquid 
nitrogen  dewar  and  heater 
used  for  temperature  con- 
•trol,  (c)  NIM  instruments 
used  for  signal  processing, 
discrimination,  counting. 


#6,218,657)  that  employs  matched  transmission  lines  to  create  an  inverted  copy  of  the 
transient  and  subtract  it  from  the  “photon+ transient”  signal,  leaving  the  photon-induced 
signal  against  a  flat,  quiet  background.  The  surrounding  electronics  must  also  include  a 
discriminator  and  amplifiers  to  generate  a  logic  output  pulse  when  the  photon-induced 
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Figure  6.  Photograph  of 
the  PCB'SPD  showing 
the  clusters  of  components 
that  implement  the  func¬ 
tions  shown  in  Figure  4. 


6 


pulse  exceeds  a  set  threshold.  In  the  IBM  prototype  QKD  system  and  in  a  separate  system 
we  subsequently  built  for  detector  characterization  [8],  these  functions  were  performed  by 
discrete  Minicircuits  RF  modules  (Figure  5(a)),  a  liquid  nitrogen  dewar  /  heater  combina¬ 
tion  for  varying  the  temperature  (Figure  5(b))  and  NIM  and  stand-alone  laboratory  instru¬ 
ments  (Figure  5(c)). 

2.  IBM  ARC  Research  Prototype  Single  Photon  Detector 

One  of  IBM’s  roles  in  this  program  has  been  to  enable  our  university  partners  to  make 
measurements  of  single-photon  detection  efficiency  on  the  novel  devices  they  are  develop¬ 
ing  that  can  be  compared  directly  with  those  made  previously  by  IBM  on  commercially- 
available  InGaAs  APDs,  so  that  the  new  devices  developed  under  this  program  can  be  fairly 
compared  with  those  used  in  earlier  work.  Transfer  to  the  universities  of  the  IBM  system  for 
single -photon  detector  characterization  proved  cumbersome  because  of  the  complexity  of 
our  original  discrete -element  laboratory  implementation.  In  order  to  remedy  this  situation, 
in  Year  Two  we  developed  an  integrated  printed-circuit-board  (PCB)  version  (Figure  6)  of 
the  single-photon  detection  system  (SPD)  that  contained  the  electronic  elements  previ¬ 
ously  implemented  using  conventional  discrete  components.  The  use  of  modem  surface - 
mount  components  and  multi-layer  PCB  techniques  made  this  a  very  compact  and  robust 
replacement  for  our  earlier  implementation. 

Although  this  system  was  originally  intended  only  for  use  by  ourselves  and  our  sub¬ 
contractors  at  UCSD  and  UTA,  great  interest  was  expressed  by  other  DARPA  QuIST  ^ 
participants  in  having  versions  of  this  system  for  their  ongoing  work  under  this  program.  In, 
response  to  this  interest,  IBM  undertook  to  provide  a  limited  number  of  research  prototype 
units  to  DARPA  QuIST  participants. 


Figure  7  (a)  PCB  attached  to  APD  module  and  fan.  The  hermetic  enclosure  surrounds  an  APD  affixed  to  a  TE  cooler,  which 
can  be  seen  in  more  detail  in  (b),  where  the  cover  of  the  module  has  been  removed 
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Figure  8.  Front  panel  of  a  dual- channel  IBM  ARC  Research  Prototype  Single  Photon  Detector. 


In  order  to  provide  a  minimum  level  of  self-sufficient  operation,  the  research  prototype 
version  of  the  PCB-SPD  included  not  only  the  PCB  itself,  but  the  APD  affixed  to  a  four-stage 
thermoelectric  cooler  and  mounted  inside  a  hermetically-sealed,  insulated  module  (Figure  7). 
One  or  two  of  these  modules  were  supplied  in  a  19-inch  rackmount  chassis  with  the  power 
supplies  necessary  to  run  the  units.  Control  and  indicators  essential  to  operation  of  the  unit 
were  brought  out  to  the  front  panel  (Figure  8).  A  rudimentary  interface  allowing  the  DC  bias 
voltage  and  the  temperature  set  point  to  be  set  remotely  via  computer  was  also  provided. 

The  completed  units  offered  a  number  of  features: 

•The  instrument  internally  provides  a  4.5-V  bias  pulse  with  length  settable  between  0.8  ns 
and  5  ns. 

•The  delay  between  the  input  trigger  and  the  bias  pulse  can  be  adjusted  over  a  range  of  —  10 
ns  in  10  ps  steps. 

•An  on-board  discriminator  produces  an  output  logic  pulse  when  the  detected  APD  signal 
exceeds  a  value-  that  can  be  set  by  the  user. 

•The  APD  can  be  cooled  to  — 210K  by  the  integrated  four-stage  thermoelectric  cooler  and 
controller. 

•A  built-in  counter  records  the  number  of  counts;  the  count  interval  can  be  set  from  1-256 
seconds  (2^  seconds,  where  0  <  n  <  8). 

•A  front  panel  meter  allows  monitoring  of  power  supply  voltages,  bias  voltage,  actual  tem¬ 
perature,  temperature  setpoint,  and  discriminator  threshold. 

•An  on-board  blanking  circuit  can  be  used  to  suppress  a  set  number  of  bias  pulses  (0-255), 
which  helps  suppress  false  counts  caused  by  afterpulsing.  (We  found  that  a  setting  of  —6, 
corresponding  to  a  12-|lIs  dead  time,  was  sufficient  to  gain  most  of  the  benefit  without 

.  unduly  increasing  the  risk  of  missing  a  desired  pulse. )  This  feature  was  added  during  a 
revision  to  the  PCB  done  during  Year  Three— it  was  not  present  in  the  basic  PCB  design 
reported  in  the  Year  Two  Progress  Report  and  shown  in  Figure  6  of  this  report. 

The  user  manual  describing  the  instrument  and  its  operation  in  greater  detail  has  been 
included  as  Appendix  A. 

We  built  a  total  of  30  single -photon  detector  modules.  Some  of  these  were  placed  into 
single  channel  instruments;  others  in  dual  channel  instruments  (such  as  that  shown  in  Fig- 
gure  8).  The  breakdown  and  disposition  of  these  units  is  given  in  Table  I. 
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The  units  delivered  to  BBN  and  Boston  University  have  already  been  integrated  into 
the  DARPA  Quantum  Cryptography  Network  in  the  Boston/Cambridge  area.  BBN  reports 
that  replacement  of  the  detectors  they  previously  used  with  the  IBM  ARC  prototype  signifi- 
cantlv  imnroved  the  nerformance  of  their  svstem  (Figure  91. 


Previous  Detectors  ARC  Detectors 


Figure  9.  Record  of  detection  events  in  the  BBN  quantum  cryptography  system.  Higher  bits  rates  with  less  noise  were  achieved 
when  the  previously-used  detectors  were  replaced  by  the  IBM  ARC  Research  Prototype  PCB-SPD.  The  recording  apparatus 
is  run  continuously;  this  explains  the  dead  time  and  variations  in  bit  rate  as  various  aspects  of  the  system  are  adjusted. 
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Figure  1 0.  Dark  counts  per  bias  puke  vs.  detection  efficiency  for  several  EPM239  APDs 


3.  Further  characterization  of  EPM  239  APDs 

The  APD  included  in  the  IBM  ARC  Research  Prototype  PCB-SPD  is  the  JDSU/Epitaxx 
EPM239.  This  detector  has  been  the  “workhorse”  most  widely  used  in  quantum  cryptography 
research  and  development  for  the  past  several  years. 

In  order  to  deliver  prototypes  with  the  best  possible  performance,  we  screened  a  large  number 
(>50)  of  EPM239’s.  The  data  gleaned  from  this  process  sheds  further  light  on  the  variability  of 
these  detectors  when  used  for  single -photon  counting  and  on  the  need  for  improvements  in  basic 
detector  technology  (as  has  been  pursued  by  UTA  and  UCSD  under  this  program). 

As  shown  in  Figure  10,  there  is  a  wide  variation  in  the  performance  characteristics  of  these 
detectors  when  cooled  to  215K  and  transiently  biased  above  breakdown.  When  used  as 
intended — at  room  temperature  and  under  appropriate  bias — the  performance  is  uniformly  good. 
This  is  not  surprising,  considering  that  these  devices  were  produced  to  meet  specifications  for 
conventional  telecom  applications  and  not  for  single -photon  counting. 

This  accumulation  of  this  much  data  on  the  performance  of  EPM239’s  enables  us  to  look  for 
correlations  between  the  performance  of  these  devices  in  single -photon- counting  mode  and  the 
parameters  that  are  traditionally  characterized  and  specified,  such  as  room-temperature  reverse 
breakdown  voltage. 

As  shown  in  Figure  11(a),  there  does  appear  to  be  a  correlation  between  dark  counts/bias 
pulse  in  single-photon  mode  and  reverse  breakdown  voltage.  According  to  Figure  1 1  (b) ,  there  is 
not  a  similar  correlation  with  afterpulse  probability,  which  seems  to  be  essentially  independent  of 
reverse  breakdown  voltage. 

The  structure  of  this  device  (Figure  12)  may  reveal  a  reason.  The  highlighted  region  shows 
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(a)  (b) 

Figure  II.  (a)  Dark  counts/bias  pulse  vs.  reverse  breakdown  voltage  corresponding  to  the  APDs  in  Figure  10.  (b) 
Afterpulse  probability  vs.  reverse  breakdown  voltage  corresponding  to  the  APDs  in  Figure  10. 

where  avalanche  multiplication  occurs.  In  the  Epitaxx  device,  the  thickness  of  this  region  is 
determined  by  a  zinc  double  diftusion  rather  than  by  epitaxy.  An  advantage  of  using  a  diffusion 
process  is  that  the  lateral  extent  of  the  region  is  also  defined,  whereas  if  the  layer  thickness  were 
defined  by  epitaxy,  the  lateral  extent  would  be  defined  by  etching  a  mesa,  which  can  create  sur¬ 
face  states  that  contribute  to  dark  current.  A  disadvantage  of  using  diffusion  is  that  this  process 
produces  variations  in  the  thickness  of  the  layer. 

In  order  to  increase  the  speed  of  devices  like  this  for  telecom  applications,  the  thickness  of  the 
layers  shown  in  Figure  12  has  been  progressively  reduced.  This  reduction  has  the  adverse  effect 
of  reducing  the  absorption,  increasing  the  probability  of  tunneling,  and  reducing  the  breakdown 
voltage,  which,  according  to  the  data  of  Figure  11,  increases  the  dark  count  rate. 
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Figure  12.  APD  device  structure  and  electric  field  profile  in  the  center  of  the  avalanche  region  (taken  from  M.  A.  Itzler, 
et  al,  "High-performance,  manufacturable  avalanche  photodiodes  for  lOGb/s  optical  receivers,"  Optical  Fiber  Communica¬ 
tions  Conference  2000,  Volume  4,  7—10  March  2000,  pp.  126-128. 


11 


4.  Conclusions 

The  IBM  ARC  Research  Prototype  PCB-SPD  has  been  provided  to  several  DARPA  partici' 
pants.  The  feedback  from  initial  use  of  these  instruments  in  quantum  cryptography  systems  has 
been  very  positive.  It  is  clear  that  the  continued  development  of  practical  quantum  cryptogra¬ 
phy  systems  depends  on  the  availability  of  good  telecom-wavelength  single-photon  detectors. 
Both  the  light-sensitive  element  itself  and  the  constellation  of  electronics  surrounding  it  must 
be  optimized  for  best  performance  in  the  single  photon  regime.  IBM  is  actively  seeking  to  license 
this  technology  to  a  partner  who  can  continue  to  make  it  available  to  the  quantum  cryptography 
community. 

B.  Separate  Hetero 'Interface  (Si/InGaAs)  Photodetectors 

The  research  of  the  UCSD  team  has  been  focused  on  four  areas  related  to  single-photon 
APD  detectors  for  quantum  cryptography; 

1.  investigation  of  theafterpulsing  effect, 

2.  development  of  physical  models  for  dark  count  and  single-photon  detection  efficiency, 

3.  demonstration  of  InGaAs-on-Si  single -photon  detectors,  and 

4.  understanding  of  APD  response  to  different  input  photon  numbers. 

Results  from  these  four  research  tasks  significantly  advance  our  understanding  of  and  ability  to 
fabricate  single-photon  APDs  at  telecom  wavelengths. 

1.  Investigation  of  afterpulsing  effect  for  single-photon  avalanche  detectors 

Summary: 

The  impact  of  the  afterpulsing  effect  and  its  temperature  dependence  were  studied  both  experimen¬ 
tally  and  theoretically.  The  primary  dark  current  and  detrap  time  as  functions  of  temperature  'were 
obtained  by  fitting  the  experimental  data.  The  result  unequivocally  reveals  the  impact  of  afterpulsing 
effect  on  the  dark  count  rate.  An  optimum  operation  temperature  for  low  dark  count  operation  was 
observed  for  single^photon  avalanche  detectors  (SPADs). 

Detailed  Discitssion; 

When  biased  above  breakdown  (Geiger  mode  condition),  an  avalanche  photodiode  (APD) 
can  produce  a  tremendous  gain  (10^-10* )  i\?ithin  a  time  window  of  a  few  nanoseconds.  Such 
Geiger-mode  APDs  must  have  a  high  detection  efficiency  for  single  photons  while  keeping  the 
dark  count  probability  low.  Unfortunately,  the  implementation  of  III-V  APDs  as  Geiger-mode 
APDs  at  wavelengths  beyond  1  /am  is  hindered  by  high  dark  count  rate.  Cooling  can  be  used 
to  reduce  the  dark  count  rate;  however,  with  InGaAs/InP  APDs  cooling  usually  worsens  the 
afterpulsing  effect.  Afterpulses  occur  when  carriers  are  released  from  traps  in  the  multiplica- 
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tion  region  and  successfully  trigger  avalanche  breakdown.  The  afterpulsing  effect  can  seriously 
degrade  the  performance  of  the  APD  and  limit  the  operation  rate  for  photon  counting. 

A  commercial  InGaAs/InP  APD  (Epitaxx  EPM239AA)  was  used  to  study  the  afterpuls- 
ing  effect.  A  laser,  synchronized  with  a  voltage  pulse  generator,  was  fired  at  a  repetition  rate  of 
50  kHz  to  generate  short  light  pulses  at  the  wavelength  of  1.3  jjm.  By  using  a  variable  attenuator, 
the  optical  energy  per  pulse  delivered  to  the  APD  was  reduced  such  that  the  average  number  of 
photons  per  pulse  was  0.3.  The  actual  number  of  photons  within  each  pulse  follows  a  Poisson 
distribution.  The  voltage  pulse  generator  sent  square  pulses  to  the  APD  at  100  kHz  (coincident 
with  and  between  the  light  pulses,  as  shown  in  Fig.  13  (a)  and  (b)).  These  pulses,  superimposed 
on  a  DC  bias,  brought  the  reverse  bias  above  Vg  for  2  ns.  The  diode  under  test  was  mounted  on 
a  temperature-controlled  plate  in  a  cryostat.  A  simple  transient- canceling  circuit  was  used  to 
separate  the  avalanche  current  pulses  from  the  capacitive  transient  pulse  at  the  diode  output. 
Finally,  the  avalanche  current  pulses  were  registered  by  a  counter. 

When  the  laser  light  was  off,  there  should  be  no  difference  between  the  count  probability  of 
coincident  and  interleaved  pulses.  The  measured  dark  count  rate  for  the  interleaved  pulses  is 
shown  in  Fig.  14  (A).  When  the  laser  light  was  on,  the  counts  for  coincident  and  interleaved 
pulses  were  registered  separately.  Since  the  coincident  pulses  had  light  input,  there  was  a  higher 
chance  of  avalanche  for  those  pulses  than  for  the  (dark)  interleaved  pulses.  The  count  rate  for 
coincident  pulses,  shown  in  Fig.  14  (♦) .  is  almost  constant  throughout  the  entire  temperature 
range  of  measurement.  What  is  most  noticeable  is  the  temperature  dependence  of  the  count 
rate  for  the  interleaved  pulses,  shown  in  Fig.  14  (□).  The  count  rate  for  the  interleaved  pulses  is 
close  to  the  count  rate  when  the  laser  light  was  off  only  at  higher  temperatures.  After  the  initial 
decreasing  trend  with  decreasing  temperature,  the  count  rate  for  the  interleaved  pulses  starts  to 
increase  when  temperature  decreases  further.  Hence,  there  is  an  optimum  operating  tempera¬ 
ture.  The  difference  between  the  dark  count  rate  and  the  count  rate  of  the  interleaved  pulses 
when  the  light  was  on  gives  a  measure  of  the  probability  of  afterpulsing  due  to  the  effect  of  carrier 
trapping. 

t _ (a)  Laser  pulses 


Coincident  pu^es  In^rleav^  pulses 

v.jyyyyyyi^ 


(b)  APD  gated  pulses 

(c)  Register  laser_on  counts 


(d)  Register  laser_on 
J  —  and  laser_off  counts 


Figure  13.  Diagram  indicating  the  different  synchronization  arrangements  used  in  the  measurements. 
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Figure  14-  Count  rates  vs.  tempera^ 
Cure.  Symbol  lines  correspond  to  mea¬ 
sured  data  using  Epitaxx  EPM239AA 
APD:  ♦  counts  registered  for  coin¬ 
cident  pulses  when  light  was  on,  □ 
counts  registered  for  interleaved  pulses 
when  light  was  on,  and  A  counts 
registered  for  interleaved  pulses  when 
light  was  off  Continuous  lines  are 
simulated  results  fitting  the  experimen¬ 
tal  data. 
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The  dark  count  probability  (Fig.  14  (A))  has  a  strong  dependence  on  the  temperature.  This 
dependence  arises  because  the  dark  counts  are  mainly  triggered  by  thermally- activated  dark  carri¬ 
ers,  and  the  activation  rate  decreases  exponentially  with  temperature.  Thus,  it  would  appear  that 
in  order  to  lower  that  component  of  the  dark  count  rate  induced  by  primary  dark  current,  the 
device  should  be  cooled  as  much  as  possible.  However,  low  temperature  exacerbates  the  after¬ 
pulsing  effect  because  the  detrap  time  of  the  carriers  increases  exponentially  with  reduced  tem¬ 
perature.  That  is  why  the  dark  count  rate  for  the  interleaved  pulses  increases  for  temperatures 
lower  than  220K  (Fig.  14  (□))•  Fitting  the  physical  model  we  developed  to  the  experimental 
data,  we  obtained  the  temperature  dependences  of  the  primary  dark  current  and  the  detrap 
time  constant  as  plotted  in  Figs.  15  and  16,  respectively. 


Figure  15.  The  dependence 
of  primary  dark  current  on 
temperature  used  in  the  curve 
fitting  shown  in  Fig.  14. 
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Figure  16.  The  dependence  of  detrap 
time  on  temperature  used  in  the 
curve  fitting  shoum  in  Fig.  14 
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2.  Development  of  physical  models  for  dark  covmt  probability  and  single-photon  quantum 
efficiency 

Summary:  q 

We  developed  a  physical  model  that  quantitatively  describes  the  behavior  of  the  dark  count  probability 
and  singk'photon  quantum  efficiency.  The  model  enables  us  to  make  direct  connections  between  the 
performance  of  SPADs  and  device  parameters  (e.g.  dark  current,  detrap  time,  gain'bandwidth  product, 
etc.)  and  operating  conditions  (e.g.  dc  gain,  Geiger  mode  gain,  gate  repetition  rate,  etc.).  Svich  infor¬ 
mation  is  useful  for  selecting  and  designing  APDs  for  single-photon  detection,  as  well  as  understanding 
device  behavior  under  operating  conditions  that  allow  single-photon  detection.  The  theory  agrees  well 
with  the  experimental  results. 

Detailed  Discussion: 

For  single-photon  detection,  a  SPAD  is  usually  dc-biased  a  few  volts  below  its  breakdown 
voltage  Fg ,  and  is  periodically  pulse -biased  above  Fg  for  a  short  time.  If  a  photoexcited  carrier 
is  present  during  this  time,  an  avalanche  process  may  be  initiated  to  produce  a  current  pulse  on 
the  order  of  milliamperes,  corresponding  to  a  Geiger  mode  gain  of  10*  -10* .  Ideally,  SPADs 
should  produce  no  current  response  under  dark  conditions.  However  in  practice,  the  seed  carri¬ 
ers  can  be  generated  not  only  by  optical  excitations  but  also  by  other  mechanisms,  such  as  pri¬ 
mary  dark  current  and  carriers  released  from  the  traps  in  the  multiplication  region.  The  primary 
dark  current  is  defined  as  the  dark  current  that  flows  through  the  same  multiplication  region  of  a 
SPAD  as  the  photocurrent  does,  measured  at  the  unity  gain  condition.  Because  dark  count  prob- 
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ability  and  single-photon  quantum  efficiency  (SPQE)  are  two  most  important  figures  of  merit 
of  SPADs,  our  analysis  will  focus  on  these  two  quantities. 

To  create  a  dark  count,  two  statistically  independent  events  have  to  happen:  (a)  dark  carri¬ 
ers  are  generated  and  (b)  the  avalanche  process  has  to  be  triggered.  Assuming  is  the  average 
number  of  dark  carriers  in  the  multiplication  region  and  P^  is  the  probability  of  triggering  the 
avalanche  event  above  the  threshold  by  each  carrier,  the  probability  of  dark  count,  P^  fol¬ 
lowing  the  Poisson  statistics  can  be  represented  as 

P,=l-exp(-Ar,/>J 


Eq.  (1)  represents  the  probability  that  at  least  one  dark  carrier  successfully  triggers  an  ava¬ 
lanche.  There  exist  four  distinct  mechanisms  in  a  SPAD  that  can  contribute  to  . 

1 )  Primary  dark  carriers  injected  or  generated  inside  the  multiplication  region  during  the  voltage  pulse. 
The  magnitude  can  be  represented  as: 

(2) 

where  is  the  primary  dark  current,  T  is  the  gate  pulse  width,  and  q  is  the  charge  of 
an  electron. 

2)  Dark  carriers  generated  before  the  voltage  pulse.  Since  the  bias  voltage  prior  to  the  pulse  is 
only  about  a  few  volts  (e.g.  1-2V)  below  ,  primary  dark  carriers  generated  before  the 
arrival  of  a  bias  pulse  can  go  through  a  series  of  impact  ionizations  with  an  average  dc  gain 
of  Mq  .  Due  to  the  finite  frequency  response  of  the  SPAD,  some  carriers  may  still  remain 
in  the  multiplication  region  when  the  voltage  pulse  arrives,  and  therefore  can  be  treated 
as  seed  carriers.  The  average  number  of  dark  carriers  thus  generated  is  given  by 


where  T*  is  the  effective  transit  time.  With  a  multiplication  gain  larger  than  the  recipro¬ 
cal  of  ionization  coefficient  ratio,  i.e.  when  M^k  >  1 ,  APDs  exhibit  a  constant  gain-band¬ 
width  product  GB.  is  then  dominated  by  the  avalanche  buildup  process,  and  can 
be  approximated  by  T*  =  /  2nGB,  where  GB  is  the  gain-bandwidth  product  of  the 

SPAD. 

3)  Dark  carriers  arising  from  the  release  of  trapped  carriers  during  the  voltage  pulse.  These  car¬ 
riers  could  have  been  trapped  during  any  of  the  previous  current  pulses  and  this  effect  is 
commonly  called  afterpulsing.  It  is  assumed  that  the  probability  of  carriers  escaping  from 
the  traps  follows  an  exponential  function  with  a  characteristic  detrap  time  constant  ; 
and  the  average  number  of  carriers  trapped  after  a  current  pulse  is  =  cM  /(l  -c)  , 
where  c  is  the  ratio  of  the  trapped  carriers  to  the  total  carriers  generated  per  avalanche 
pulse  and  its  value  depends  on  the  number  of  unfilled  traps  in  the  multiplication  region. 
The  average  number  of  detrapped  carriers  within  the  duration  of  the  voltage  pulse  t  is 
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where  AT  is  the  reciprocal  of  the  pulse  repetition  rate. 

4)  Based  on  a  argument  similar  to  that  in  2),  carriers  released  from  the  traps  before  the  voltage 
pulse  may  also  stay  in  the  multiplication  region  when  the  voltage  pulse  arrives.  The  number  of 
such  dark  carriers  can  be  written  as 


exp 


^  nAT 


-exp 


^  nAT^ 


=  P,K 


exp(r;/T^)-l 
exp(A7' /Tj)-1 


(5) 


Finally,  the  total  number  of  dark  carriers  per  pulse  in  a  SPAD  is  the  sum  of  these  four  sources,  so 


^cl  ~  ^DMl  ^DM2  -^(1  ^12 


(6) 


Another  important  parameter  for  the  dark  count  calculation  is  the  probability  for  a  carrier 
to  cause  an  avalanche,  P^.  McIntyre  has  theoretically  calculated  P^  ’s  of  electrons  and  holes  as 
functions  of  the  ionization  coefficient  ratio  k  and  5  =  ^a  dx  {a  :  ionization  coefficient  for  elec¬ 
trons).  Here  we  treat  as  a  parameter  so  our  results  are  applicable  to  any  APDs  without  being 
restricted  to  a  particular  device  structure  or  electric  field  distribution  in  the  multiplication  region. 

Substituting  Eqs.  (2)-(6)  into  Eq.  (1),  the  dark  count  probability  P^  is  given  by 


P, 


1-exp  •{ 


-P 


exp(T/T,)-l  exp(r-„fa,)-l 

q  InqGB  ‘‘  exp(Ar/T^)-l  exp(A7’/Tj)-l 


(7) 


Because  Pj  appears  in  both  sides  of  Eq.(7),  this  equation  must  be  solved  numerically. 


The  SPQE,  also  known  as  detection  efficiency,  is  commonly  defined  as 

SPQE  =  {P„„-P,)IP^„  (8) 

where  is  the  probability  of  a  current  pulse  triggered  by  either  a  photon-generated  carrier  or 
a  dark  carrier  when  the  single-photon  source  is  on;  and  =  1  -  e~^°  is  the  probability  that 
an  incident  optical  pulse  contains  at  least  one  photon  with  the  average  number  of  photons  per 
pulse  Nq  .  Although  this  is  a  commonly  accepted  definition  of  SPQE,  we  should  point  out  that 
Eq.  (8)  describes  the  true  single -photon  quantum  efficiency  only  when  the  error  rate  due  to  after¬ 
pulse  is  relatively  low.  Otherwise,  Eq.  (8)  tends  to  overestimate  SPQE. 

The  method  of  calculating  P^„  is  similar  to  the  calculation  of  the  dark  count  probability,  P^ , 
except  that  the  effect  of  optical  input  has  to  be  considered.  Accordingly,  Eqs.  (1)  and  (6)  should 
be  modified  as 


Figure  j 
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(9) 


N  -N  +N  +N'  -^N'  +nN 

DM\  ^  DM2  ^  t\  ^  il 


(10) 


where  rj  is  the  conventional  quantum  efficiency  of  a  detector,  determined  by  the  absorption 
coefficient  and  the  thickness  of  the  absorption  layer.  Equations  (4)  and  (5)  can  still  be  used  for 
N'^  and  N'2  except  that  Pj  should  be  replaced  by  The  additional  term  in  Eq.  (10), 
accounts  for  the  optical  induced  carriers.  The  final  expression  of  is  then  given  by 


1-exp 


-P 


exp(.fa)-l 

on  it  _ /  A  'T’  / _  \  1  " 
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exp(A7’/T^)-l 


exp(TV/T^)-l 

exp(Ar/Tj-l 


+  riK) 


(11) 


Finally,  the  SPQE  can  be  calculated  from  Eqs.  (8)  and  (11). 

The  above  analysis  shows  explicitly  how  dark  count  and  SPQE  depend  on  different  device 
parameters  and  operating  conditions.  In  the  following,  we  will  use  calculated  results  to  illustrate 
the  effects  of  the  major  parameters.  The  dependence  of  Pj  as  a  function  of  SPQE  is  shown  in 
Figure  17  for  different  (0.1,  1  and  10  pA),  assuming  that  afterpulsing  is  negligible,  which  is 
reasonable  given  «:  AT  .  Note  the  strong  effect  of  primary  dark  current  on  the  dark  count 
probability.  The  calculation  agrees  very  well  with  the  measured  results  obtained  from  commercial 
APDs  (EG&.G  InGaAs  APD  at  240  K  with  measured  of  5pA  and  NEC  NDL3351P  APD 
at  130  K  with  measured  of  0.1  pA),  as  shown  in  Fig.  17.  The  same  measurement  condi¬ 

tions  were  used  in  the  calculation  shown  in  Fig.  17.  In  the  region  of  high  quantum  efficiency,  the 


Figure  1 7.  Calculated  (solid 
lines)  and  measured  (symbol 
lines)  dark  count  probability 
vs.  single^photon  quantum  effi¬ 
ciency  for  different  Other 
parameters  are:  /  =  1  /  AT 
=  I  OOkHz,  pulse  width  t  = 
2ns,  '^d  =  200  ns,  cis  1%, 

Ng  =  0.3  photon/pulse,  M  = 
10^,andGB  =30GHz. 
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measured  dark  count  increases  more  rapidly  than  the  model.  This  is  because  the  dc  bias  voltage 
and  therefore  increases  in  the  measurement  in  order  to  reach  the  high  SPQE,  whereas  our 
calculation  has  kept  the  dc  bias  and  Mg  constant. 

The  effect  of  the  dc  gain  is  shown  in  Fig.  18  for  SPADs  with  different  gain-bandwidth 
products.  Here  the  is  assumed  to  be  O.lpA,  and  all  the  other  parameters  are  kept  the  same 
as  in  Fig.  17.  Mg  affects  through  as  shown  in  Eq.  (3)  -  a  higher  Mg  yields  a  larger  P^ . 

The  results  from  Figs.  17  and  18  suggest  that  low  primary  dark  current  and  low  dc  bias  gain  are 
important  for  achieving  a  low  dark  count  or  equivalently  a  high  single-photon  quantum  efficiency  ‘ 
at  a  given  dark  count. 


Figure  1 8.  Calculated 
dark  count  probability 
vs.  single-photon  (Quan¬ 
tum  efficiency  for  three 
different  Mg  levels.  The 
solid  lines  indicate  APDs 
with  GB  of  30  GHz; 
the  dashed  lines  indicate 
APDs  with  GB  of  100 
GHz-  The  detrap  time 
constant  is  set  much 
smaller  than  the  time 
period  between  pulses,  so 
that  the  afterpulse  effect 
is  almost  negligible  in  this 
calculation. 


Single-photon  quantum  efficiency  % 


Experimental  results  have  also  shown  that  the  afterpulse  contribution  to  the  measured  counts 
is  negligible  at  low  pulse  repetition  rate.  This  can  be  clearly  explained  from  the  last  two  terms  in 
the  product  of  Eq.  (7).  When  all  the  trapped  carriers  are  released  before  the  next  voltage  pulse 
arrives,  the  chance  of  afterpulsing  is  minimal.  Experimentally,  Si  SPADs  have  exhibited  weak 
trapping  with  very  small  (~  10  ns  at  room  temperature) .  In  contrast,  for  IIW  APDs  is  typi¬ 
cally  hundreds  of  nanoseconds  or  microseconds  at  room  temperature  and,  in  addition,  the  latter 
devices  have  a  higher  trap  density  in  the  IIIA7  multiplication  region.  As  the  device  is  cooled  to 
achieve  low  primary  dark  current,  for  IIW  APDs  increases,  causing  a  more  serious  afterpuls¬ 
ing  effect.  In  Figure  19(a),  the  relative  contribution  of  afterpulse  to  the  total  dark  count  is  plot¬ 
ted  against  the  pulse  repetition  rate  with  different  detrap  times,  .  Clearly  the  contribution 
of  afterpulse  increases  rapidly  and  becomes  dominant  once  the  repetition  rate  reaches  a  certain 
threshold  value.  Figurel9(b)  plots  the  voltage  pulse  period  (=  1 //)  for  which  the  afterpulse 
contribution  reaches  1%  of  the  total  dark  count  versus  for  different  M^ .  We  conclude  that  to 
keep  the  afterpulse  at  a  negligible  level,  two  conditions  should  be  satisfied:  (1)  the  voltage  pulse 
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(a) 


(b) 


Figure  19.  (a)  Contribution  of  afterpubes  to  dark  counts  vs.  pulse  repetition  rate  for  different  detrap  time  constants  with 

=  l<y.  (b)  Voltage  pube  period  (Iff)  at  which  the  contribution  of  afterpubes  to  dark  counts  is  1%  as  a  func¬ 
tion  of  ,  with  different  Geiger  mode  gains.  =  O.lpA,  SPQE=20%  and  all  the  other  parameters  are  the 
same  as  Fig.  17.. 
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period  should  be  at  least  ten  times  longer  than  the  detrap  time  tj ;  (2)  a  low  pulse  repetition  rate 
is  needed  when  a  higher  Geiger  mode  gain  is  used  for  single -photon  detection. 

3.  High-performance  InGaAs-on-Si  single  photon  avalanche  photodetectors 
Summary: 

We  demonstrated  an  InGaAs'OmSi  single  photon  avalanche  diode  (SPAD)  with  high'Singk'photon 
quantum  efficiency  and  low- dark' count  probability  under  gated'mode  operation.  At  223  K,  the  device 
exhibited  a  dark  count  probability  o/ 10"*  -  10'^  and  a  singC'photon  quantum  efficiency  of  33%  without 
antUreflection  coating.  The  devices  also  showed  high  efficiency  and  low  dark  count  with  the  smallest 
normalized  excess  voltage  (4-  7%  ,  as  a  result  of  the  Si  multiplication  region 

Detailed  Discussion: 

Figure  20  schematically  depicts  the  Separated  Absorption  and  Multiplication  (SAM)  struc¬ 
ture  of  the  InGaAs-on-Si  APDs.  The  device  has  a  light-sensitive  diameter  of  150  ju.m  and  a  total- 
junction  diameter  of  200  /am. 

hv  ) 

(  p-metal  contact  pad 


n-metal  contact  pad 

n-Si  contact  layer 


Figure  20.  Schematic 
InGoAs-ori'Si  SPAD  device 


structure. 


p-!nP  contact  layer 


InGaAs  (absorption) 


Si  (multiplication) 


Figure  2 1  shows  typical  1-V  characteristics  in  the  dark  and  under  illumination  at  298  K  and 
223  K.  At  223K  the  dark  current  at  90%  of  the  breakdown  voltage  is  as  low  as  20  pA.  For  the 
single  photon  counting  experiment  the  devices  were  connected  on-wafer  with  high-frequency 
probes  and  operated  in  a  gated  mode  where  a  train  of  voltage  pulses  was  superimposed  on  a  DC 
bias,  The  DC  bias  was  set  slightly  below  the  device  breakdown  voltage.  The  voltage  pulse 
was  2  ns  wide  and  1.9  V  high  and  had  a  repetition  rate  of  100  kHz.  A  190  ps  optical  pulse  was 
generated  using  a  laser  pulse  source  (Hamamatsu  C8898)  emitting  at  a  wavelength  of  1550  nm. 
The  laser  pulses  were  attenuated  to  an  average  energy  of  0.1  photon/pulse  and  coupled  to  the 
device  with  a  single -mode  optical  fiber  probe.  The  avalanche  pulses  produced  by  the  SPAD, 
caused  by  photo  or  dark  carriers,  were  registered  with  an  electronic  counter  (Stanford  Research 
Systems  SR400) .  A  simple  voltage  transient-canceling  circuit  was  used  to  separate  the  avalanche 
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current  pulses  from  the  capacitive  transient  pulses  at  the  diode  output.  Both  the  pulsed  laser  and 
the  counter  were  triggered  by  the  voltage  pulse  generator  (Agilent  81110A).  The  superimposed 
gated  bias,  +  F^  ,  produced  an  excess  voltage  F^  above  the  APD  breakdown  voltage,  F^,^ ,  fol¬ 
lowing  the  simple  relation:  ~  Ka  •  experiment,  the  height  of  the  voltage  pulse 

Fg  remained  fixed  at  1.9  V  and  the  excess  voltage  was  varied  by  changing  the  DC  bias.  While 
was  varied  from  -37.6V  to  -38.1V,  the  excess  voltage,  ,  was  changed  from  3.4%  to  4.7%  F^^ 
(F,^= -38.2  V  at223K). 


Figure  2 1 .  Photo  and  dark  current  versus 
reverse  bias  of  InGoAs-on-Si  SPAD  at  dif¬ 
ferent  temperatures. 


Dark  count  probability  ( )  versus  single-photon  quantum  detection  efficiency  {SPQE)  of 
the  InGaAs-on-Si  SPAD  is  plotted  in  Figure  22.  Here  P^  is  defined  as  the  probability  of  having 
a  count  triggered  by  non-photon- excited  carriers  during  a  gate  period.  P^  was  obtained  from 
the  ratio  of  the  number  of  registered  counts  and  the  total  number  of  gate  pulses  when  the  device 
was  in  the  dark.  The  SPQE  is  defined  as  the  probability  of  having  a  count  triggered  by  photo¬ 
generated  carriers  during  a  gating  period.  At  each  excess  bias  level,  the  SPQE  was  obtained  by  1) 
subtracting  the  dark  count  from  the  total  number  of  registered  counts  when  the  laser  was  turned 
on  and  2)  normalizing  the  result  by  the  number  of  available  photons  per  pulse.  Because  the  light 
source  was  a  semiconductor  laser,  a  Poisson  distribution  for  photon  number  was  used  to  calcu¬ 
late  the  average  number  of  photons  per  light  pulse.  In  our  measurements  the  afterpulse  effect 
was  negligible  because  the  gate  repetition  rate  (100  kHz)  was  chosen  to  be  much  lower  than  the 
escape  rate  of  trapped  carriers.  If  the  gate  repetition  rate  was  increased  substantially,  the  above 
method  would  overestimate  the  SPQE  due  to  the  afterpulse  effect. 


Each  data  point  in  Fig.  22  was  measured  at  a  corresponding  DC  bias.  Both  P^  and  SPQE 
increase  with  the  magnitude  of  the  DC  bias.  When  the  DC  bias  voltage,  F^^ ,  reaches  the  break- 
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down  voltage  the  dark  count  probability,  Pj ,  experiences  an  abrupt  jump,  limiting  the  maximal 
SPQE  that  the  device  can  practically  achieve.  Fig.  22  shows  that  the  dark  count  probability 
was  in  the  range  of  8  x  10'^  -  4x  10'^  when  the  SPQE  was  between  2.5%  and  33%,  for  a  repeti¬ 
tion  rate  of  100  kHz  at  223K.  This  result  is  comparable  to  the  published  results  on  InGaAs/InP 
APDs  under  similar  conditions.  Note,  however,  that  the  commercial  APDs  used  for  comparison 
had  a  smaller  diameter  (40  ixm) .  Small  APDs  are  expected  to  exhibit  low  dark  count  probabili¬ 
ties  at  the  same  SPQE  .  Because  the  area  of  the  InGaAs-on-Si  SPAD  is  25  times  larger  than  of 
the  commercial  InGaAs/InP  APDs,  it  is  expected  to  have  a  more  than  one  order  of  magnitude 
improvement  if  we  decrease  the  device  size  down  to  40  ptm.  Also,  a  40%  improvement  in  SPQE 
could  be  achieved  using  an  antireflection  coating.  Finally  it  is  possible  that  the  results  may  have 
been  negatively  affected  by  the  fact  that  the  measurements  were  performed  on-wafer  instead  on  a 
packaged  device. 

Values  for  the  noise  equivalent  power  (NEP)  based  on  single -photon  detection  were  also 
obtained.  Its  dependence  on  the  discrimination  levels  of  the  pulse  counter  is  shown  in  Fig.  23.  It 
appears  that  the  NEP  remains  nearly  constant  with  the  DC  bias,  ,  up  to  a  certain  discrimination 
level.  Above  this  level  the  NEP  increases  dramatically.  The  minimum  NEP  for  different  values 
of  falls  in  the  range  of  8-18x  10'‘*  W/Hz'^,  which  is  comparable  to  the  published  values  for 
InGaAs/InP  SPADs  under  similar  conditions. 

One  noteworthy  feature  of  the  InGaAs-on-Si  Geiger-mode  SPADs  is  that  the  required  excess 
bias  to  achieve  the  highest  detection  efficiency  of  33%  was  only  1.9V  (=4.7%  To  the  best 
of  our  knowledge  this  is  by  far  the  smallest  excess  bias  required  to  obtain  such  high  SPQE  while 
keeping  the  dark  count  low.  The  ability  to  simultaneously  achieve  high  SPQE  and  low  dark  count 
probability  at  low  excess  bias  is  highly  desirable  for  practical  systems,  because  it  relaxes  the  require¬ 
ments  for  bias  stability  and  temperature  control,  as  well  as  improves  the  device  reliability.  This 
resultd  agrees  with  theoretical  analysis,  which  indicates  that  APDs  with  small  k  values  have  higher 
avalanche  breakdown  probability  under  a  given  normalized  excess  bias,  . 
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Figure  23.  Noise  equivalent  power 
versus  discrimination  level  for 
InGoAs-on-Si  SPAD  at  different 
DC  biases. 


Discrimination  level  (mV) 


4.  APD  response  to  dif¥erent  numbers  of  input  photons 

Summary: 

For  the  first  time  we  obtained  the  response  of  Geiger-mode  APDs  for  different  numbers  of 
photon  or  dark  carriers.  We  have  shown  the  avalanche  probabilities  for  one,  two,  and  three  pho¬ 
tons  at  different  discrimination  levels. 

Detailed  Discussion: 

Although  it  has  been  shown  that  avalanche  photodiodes  (APDs)  working  in  Geiger  mode  are 
able  to  detect  single  photons,  there  have  been  few  studies  on  how  Geiger-mode  APDs  respond  to 
different  numbers  of  input  photons.  The  lack  of  such  studies  has  been  due  partly  to  the  notion 
that  the  APD  output  might  be  too  noisy  to  render  any  useful  information  on  photon  number  and 
partly  to  the  unavailability  of  light  sources  able  to  generate  specific  photon  numbers  accurately. 
We  have  developed  a  method  that  can  yield  avalanche  probabilities  for  Geiger-mode  APDs  under 
different  numbers  of  photon  excitation  using  the  dark  count  measurement  and  a  mathematical 
routine. 

We  assumed  that  the  APD  response  to  a  dark  carrier  is  equivalent  to  its  response  to  a  photo¬ 
carrier  generated  at  the  same  time.  Furthermore,  since  the  dark  count  probability  was  typically 
much  lower  than  1%,  we  assumed  that  all  the  dark  counts  were  caused  by  single  carrier  and 
therefore  the  APD  response  to  dark  carriers  can  represent  its  response  to  single  photo-carrier. 

The  APD  is  biased  below  its  breakdown  voltage  and  pulsed  above  breakdown  for  a  few  nanosec¬ 
onds,  under  the  so-called  gated  Geiger  mode  operation.  The  data  presented  here  were  obtained 
from  a  commercial  InGaAs/InP  APD  that  was  biased  below  breakdown  and  pulsed  above  break¬ 
down  for  2  ns  with  various  pulse  voltages  at  a  rate  of  lOOkHz.  The  output  avalanche  coming  out 
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of  the  APD  was  fed  into  a  pulse  counter.  By  measuring  the  count  number  at  different  discrimina¬ 
tion  voltage  levels,  F^,.^ ,  the  output  distribution  of  the  Geiger-mode  APD  was  obtained. 

Figure  24  shows  the  measured  dark  count  distribution  at  room  temperature  against  the  dis¬ 
crimination  voltage.  The  dark  count  probability  can  be  expressed  as  P^  =  \- exp ) ,  where 
is  the  number  of  dark  carriers  within  the  2  ns  pulse  that  may  trigger  avalanche  events,  and  P^^ 
is  the  single-carrier  initiated  avalanche  probability  using  as  the  measurement  criteria.  Since 
is  given  by  /  q ,  where  7^^  and  t  are  the  primary  dark  current  and  pulse  width,  one  can 
directly  calculate  the  “single-carrier  induced  avalanche  probability”  7^,  =-9ln(l-7^)/(7^^T)as  a 
function  of  from  the  data  in  Fig.  12.  The  probability  density  function  (pdf)  of  a  single -carrier 
initiated  avalanche,  ,  is  obtained  by  differentiating  7^,  with  respect  to  .  One  major  step 
here  is  to  derive  the  avalanche  probability  and  pdf  for  2-carriers,  3-carriers,  {dP^^^dP^^  )  and  so  on 
as  functions  of  .  Because  every  carrier  contributing  to  the  output  current  must  be  the  offspring 
of  one  of  the  carriers  that  initiate  the  avalanche,  we  can  derive  the  following  relation  (Eq.  12). 

‘‘P,  =  1"  dP.,  (vW,..„(F  -  V)  dvj  -  v)  rfv  rfl' 

From  Eq.  12,  we  can  further  calculate  <77^„dPan  as  long  as  we  know  all  dP^/s  where 
/■  =  1, 2, . . . , «  - 1 .  Since  we  have  obtained  7^,  and  dP^^  .from  the  dark  count  measurement  as  stated 
before,  we  can  hereby  obtain  Geiger-mode  APD  response  to  any  number  of  input  carriers.  The 
avalanche  probabilities  density  functions  and  avalanche  probabilities  caused  by  1,  2,  and  3  carriers 
are  plotted  against  in  Figures  13  and  14. 

To  our  best  knowledge,  this  is  the  first  time  that  the  relation  between  avalanche  probability 
distribution  and  photon  number  is  derived.  For  practical  apg^cations,  this  study  shows  us  how  to 
trade  off  APD  sensitivity  against  dark  count  by  setting  the  discrimination  voltage.  In  the  particular 


Figure  24.  Measured  dark 
count  probability  as  a  function 
of  discrimination  level  at  room 
temperature.  is  the  excess 
pulse  height  on  top  of  APD’ s 
breakdoum  and  is  the  APD 
breakdoum  voltage. 
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case  of  Fig.  26,  for  example,  one  can  set  at  70  mV  so  the  Geiger-mode  APD  has  about  18% 
avalanche  probability  for  3  or  more  photocarriers  but  negligible  dark  counts  since  the  probability  of 
generating  more  than  one  dark  carrier  within  the  narrow  time  window  even  at  room  temperature 
is  very  low. 


Figure  25.  Probability  density  (dP ) 
of  one-  (solid  lines),  tuio-  (dash  lines) 
and  three-(dot  lines)  carrier  initiated 
avalanche  versus  discrimination  level, 
derived  from  the  dark  count  measure¬ 
ment  (Fig.  24)  and  Equation  (12). 


Figure  26.  Dependence  of  avalanche 
probability  on  discrimination  levels 
for  one- (solid  lines) ,  two-  (dash 
lines)  and  three  (dot  lines) -carrier 
initiated  avalanche. 
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C.  Optimization  of  low-noise  InGaAs  APDs  for  detection  of  single  photons  at  telecom 
wavelengths  (UT  Austin) 

In  this  reporting  period  we  have  focused  on  a  study  of  the  origin  of  dark  counts  in  an 
In^  gjGag^^  As  /  In^  ^gAs  SACM  APD. 

The  type  of  photon- counting  APD  studied  was  a  SACM  (Separate -Absorption-Charge  Mul¬ 
tiplication)  APD  with  a  1  jU,m- thick  In^ggCa^^^As  absorbing  layer  and  400nm- thick  In^  g^Alg  ^gAs 
multiplication  layer.  A  schematic  of  the  detector  is  shown  in  Figure  27.  The  structure  was  grown 
on  n"*"  substrate  using  metal-organic  chemical  vapor  phase  deposition  (MOCVD).  Top-illuminated 
mesa  devices  were  fabricated  using  standard  photolithography,  lift  off,  wet  chemical  etching  and 
SiOj  passivation.  Ti/Pt/Au  and  AuCe/Ni/Au  were  used  as  the  p-contact  and  n-contact,  respec¬ 
tively.  Details  of  the  structure  and  the  processing  can  be  found  in  Ref.  [1]. 
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Figure  27.  Schematic  of  /  In^j^Alp^gAs  APD 

The  device  was  investigated  using  the  gated  mode  of  operation.  Figure  28  shows  an  Arrhenius 
plot  of  dark  count  rate  at  different  excess  voltages.  These  values  of  excess  bias  yield  single  photon 
detection  efficiencies  in  the  range  ~  10-15%.  Since  the  dark  counts  exhibit  an  exponential 
dependence  on  temperature,  an  activation  energy  can  be  extracted  from  the  slope  of  a  plot  of  the 
logarithm  of  the  dark  count  rate  versus  1/  kT .  The  activation  energies  obtained  from  the  plots 
range  from  0.12  to  0.15eV  for  excess  bias  of  3.5V  to  1.5V. 

The  value  of  activation  energy  can  be  used  to  identify  the  responsible  mechanism.  The  fol¬ 
lowing  mechanisms  are  possible  thermal  sources  of  dark  counts  in  the  SACM  APD: 

1)  Dark  counts  resulting  from  diffusion  of  carriers  from  quasi-neutral  n-  and  p- 
regions  to  the  multiplication  layer  have  an  exponential  dependence  on  tem¬ 
perature  through  the  band  gap.  Flence  the  activation  energy  is  expected  to  be 
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equal  to  the  band  gap  of  the  neutral  region. 

2)  For  dark  counts  resulting  from  generation-recombination  in  the  In^  gjOa^  As 
absorber  layer,  the  activation  energy  is  expected  to  be  ~  Vi  the  band  gap  of 
Im  „Ga,,,As,  which  is  ~  0.35eV. 

3)  For  generation-recombination  in  the  multiplication  layer,  the  activation  energy 
should  be  ~  Vi  band  gap  of  In^  j^Alg^As,  which  is  0.6 7 eV. 

The  experimentally- observed  activation  energy  is  much  lower  than  the  value  expected  from 
any  of  the  mechanisms  mentioned  above.  In  the  following,  we  explain  the  temperature  depen¬ 
dence  of  dark  count  rate  by  invoking  hand-to-hand  tunneling  in  the  multiplication  layer.  Band 
profile  simulations  of  the  device  above  breakdown  voltage  indicate  a  triangular  barrier  for  tunnel¬ 
ing. 


Figure  28.  Dark  count  rate  as  a  function  of  1/kT  for  different  excess  voltages. 

Dark  counts  due  to  tunneling  ( )  can  be  estimated  by  dividing  the  tunneling  dark 
current,  I,umeimg  i  by  the  electronic  charge,  q,  and  multiplying  by  the  breakdown  probability, 
PehiJunction) . 

^d,ur,„di„g  =  P Amotion) 

(1) 

PAJunction)  is  the  electron-hole  pair  breakdown  initiation  probability  at  the  junction, 
where  the  electric  field  is  maximum[2].  A  more  thorough  analysis  involves  integrating  the  prod 
uct  of  the  generation  rate  due  to  band-to-band  tunneling  and  position- dependent  breakdown 
probability  over  the  depletion  layer  width. 
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The  dark  counts  per  second  due  to  tunneling  can  be  written  as  [3] : 


N. 


d,fmmeling 


=  Peh  {junction)  ■  exp 


-a 


C,  can  be  expressed  as 


C.=, 


|2- ffj. 


A-n^-h^ 


(2) 


(3) 


The  form  of  C,  depends  on  the  type  of  tunneling  barrier.  For  a  triangular  barrier, 


Cj  =  exp 


-3/2  ^ 
g 


3-q-h-E„, 


(4) 


The  various  terms  in  equations  (2) -(4)  are  described  below.  E^  is  the  electric  field  in  the 
multiplication  region.  Equation  (2)  is  derived  assuming  a  constant  electric  field  in  the  depletion 
layer  [4] .  Since  the  electric  field  in  the  multiplication  layer  of  an  SACM  APD  is  not  constant,  the 
average  electric  field,  obtained  by  solving  Poisson’s  equation  for  the  SACM  APD  with  appropriate 
boundary  conditions,  was  used  in  the  calculation.  is  the  voltage  drop  across  the  junction,  E^ 
is  the  band  gap  of  In^jjAlg^gAs,  h  is  Planck’s  constant,  and  w*  is  the  tunneling  effective  mass.  Its 
value  is  taken  to  be  the  reduced  effective  mass  defined  as  follows, 

I  1  1 


m. 


m. 


m,. 


(5) 


where,  w*  =  0.079  is  the  T  valley  electron  effective  mass  and  =  0.086  is  the  light  hole 
effective  mass  for  In^j^Alg^gAs.  Note  that  equation  (2)  applies  only  to  direct  tunneling,  which  is 
dominant  for  direct  band  gap  In^g^Alg^gAs. 


Figure  28.  Experimental  and  calculated  dark  count  rate  as  a  function  of  temperature  for  excess 
bias  of  1.5V. 
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Figure  28  shows  the  calculated  dark  count  rate  as  a  function  of  temperature,  normalized  to  the 
value  at  room  temperature.  It  should  be  noted  that  the  absolute  value  of  dark  count  rate  depends 
on  experimental  details:  for  example,  the  discriminator  threshold,  the  quenching  technique,  and 
the  repetition  rate.  Also,  the  simple  tunneling  model  does  not  take  into  account  non-uniformities 
in  the  active  area.  However,  regardless  of  the  absolute  value  of  the  dark  counts,  for  a  given  excess 
voltage  their  temperature  dependence  is  unique.  The  temperature  dependence  of  dark  count  rate 
is  primarily  due  to  the  temperature  dependence  of  the  band  gap  of  In^jjAlg^gAs.  Figure  29  shows  a 
plot  of  the  of  lUg  gjAlg  ^gAs  versus  MkT ,  using  the  temperature  dependence  of  the  band  gaps 
of  InAs  and  AlAs  and  the  bowing  parameter  for  In^g^Alg^As.  In  the  temperature  range  of  200K  to 
300K,  has  an  almost  linear  dependence  on  M kT  with  the  slope  ct  =  4.8x10  ^  eV^'^ .  Using 
E^J^  as  a  linear  function  oi  \lkT  plot  of  the  log  of  dark  count  rate  versus  MkT  has  a  slope  of 
-a 

V  / 

The  activation  energy  expected  from  this  expression  at  an  excess  bias  of  3.5V  is  0.14eV,  which 
matches  reasonably  well  with  the  experimentally  observed  value  of  0.1  leV.  The  activation  energy 
depends  inversely  on  electric  field  in  the  multiplication  layer,  which  was  also  observed  experimen¬ 
tally. 


1/kT(1/eV) 

Figure  29.  Dependence  of  on  I  /kT 


It  is  interesting  to  note  that  in  this  device,  even  though  the  dark  count  generation  mecha¬ 
nism  is  dominated  by  tunneling,  it  is  not  the  primary  dark  current  generation  mechanism  below 
breakdown.  Figure  30  shows  the  dark  current  below  breakdown  voltage  for  the  device  at  various 
temperatures.  For  a  broad  range  of  voltages  above  the  punch  through  value,  the  activation  energy 
is  0.4  eV.  This  points  to  generation- recombination  in  In^ggOag^^As  as  the  responsible  mechanism. 
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Reverse  bias  (V) 

Figure  30.  Dark  current  as  a  function  temperature 

The  dominance  of  hand-to-hand  tunneling  in  Geiger  mode  operation  for  this  device  is  attrib¬ 
uted  to  high  electric  field  in  its  thin  multiplication  layer  of  400  nm.  However  such  thin  multipli¬ 
cation  layer  results  in  very  good  timing  performance.  A  75-/.tm  diameter  device  exhibits  timing 
jitter  of  only  60  ps  at  200  K  when  operated  at  3.5V  above  breakdown  voltage.  Thus  there  is,  a 
trade  off  between  the  timing  performance  of  a  Geiger  mode  APD  and  tunneling  in  its  multiplica¬ 
tion  layer.  Since  tunneling  current  is  very  sensitive  to  electric  field,  we  estimate  that  increasing 
the  multiplication  layer  thickness  by  a  factor  of  2  will  reduce  the  tunneling  induced  dark  counts 
by  3  orders  of  magnitude.  The  increased  thickness  of  the  multiplication  layer  is  also  expected  to 
give  a  steeper  increase  in  breakdown  probability  with  excess  voltage,  and  hence  improved  single 
photon  counting  performance  [5].  Ultimately  generation-recombination  in  the  multiplication 
layer  and  avalanche  build  up  time  will  determine  the  optimum  multiplication  layer  thickness. 

In  summary,  the  origin  of  dark  counts  in  an  In^  jjGa^  As  /  In^  Alp  ^gAs  S  ACM  APD  was 
studied  using  the  temperature  dependence  of  the  dark  count  rate.  The  dark  count  rate  in  these 
APDs  was  found  to  be  dominated  by  band-to-band  tunneling  in  the  multiplication  layer. 

D.  Detection  of  Single  Photons  Using  Sum  -Frequency  Generation 

Another  approach  to  detection  of  single  photons  at  telecom  wavelengths  that  has  been 

pursued  under  this  program  is  to  use  sum-frequency  generation  to  increase  the  energy  of 

the  single  photon  so  that  a  silicon  APIKeHyEbgCu<(dd2Haerigt6lftii6'n.  The  proposed  imple- 

mention  uses  a  periodically-poled  W#ti^uide  to  produce  efficient 

sum-frequency  mixing  of  a  strong  a  weak  C-band  signal 
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(—1550  nm)  to  produce  a  signal  at  the  sum  frequency  (—710  nm),  which  can  be  detected 
using  a  Geiger-mode  silicon  APD.  Overall,  the  effort  devoted  to  this  approach  has  been  less 
than  originally  planned,  because  as  the  QulST  program  evolved,  we  believed  that  pursuing 
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Figure  3 1 .  Experimental  set-up  for  testing  the  sum-frequency  generation  approach  to  single  photon 
detection. 

the  integrated  detector  electronics  described  above  represented  a  more  immediate  contribu¬ 
tion  to  this  program  and  to  the  field  of  quantum  key  distribution.  However,  in  earlier  years 
we  made  some  preliminary  measurements  (such  as  the  sensitivity  of  silicon  detectors  at 
telecom  wavelengths) .  In  Year  Three  of  the  program,  we  assembled  and  tested  a  preliminary 
set-up  for  investigating  sum-frequency  generation  (Figure  31). 


Wavelength  (nm) 


Figure  32.  Scans  of  sum-frequency  output  versus  O-band  wavelength. 


The  arrangement  shown  in  Figure  31  allowed  light  in  both  the  O-band  and  C-band  wave 
length  ranges  to  be  coupled  simultaneously  into  the  PPLN  waveguide.  A  fiber  with  a  lensed 
end  was  used  for  efficient  coupling  into  the  waveguide,  which  was  mounted  on  a  tempera¬ 
ture-controlled  stage.  Both  the  temperature  of  the  PPLN  crystal  and  the  wavelength  of  the 
input  light  could  be  tuned  in  order  to  achieve  phasematching. 

Figure  32  shows  the  variation  of  the  intensity  at  the  sum-frequency  as  the  wavelength  of 
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the  O'band  input  is  scanned  while  the  wavelength  of  the  C-band  input  is  held  constant  at 
1550  nm.  Phasematching  for  sum-frequency  generation  occurs  when  the  shorter-wavelength 
is  tuned  to  1317.6  nm,  producing  light  at  a  sum-frequency  of  712  nm. 

Preliminary  characterization  of  these  waveguides  revealed  several  unanticipated  aspects 
of  their  behavior.  Although  a  sum-frequency  signal  was  detected,  self-doubling  of  the  O- 
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band  signal  was  considerably  stronger.  This  behavior  was  unexpected,  given  the  design  of  the 
PPLN  waveguide,  which  should  not  provide  phasematching  for  a  strong  self-doubling  inter¬ 
action.  We  are  continuing  to  investigate  the  cause  of  this  behavior. 
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